Non-monotonic density dependence of the diffusion of DNA fragments in low-salt 

suspensions 
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The high linear charge density of 20-base-pair oligomers of DNA is shown to lead to a striking non- 
monotonic dependence of the long-time self-diffusion on the concentration of the DNA in low-salt 
conditions. This generic non-monotonic behavior results from both the strong coupling between 
the electrostatic and solvent-mediated hydrodynamic interactions, and from the renormalization 
of these electrostatic interactions at large separations, and specifically from the dominance of the 
far-field hydrodynamic interactions caused by the strong repulsion between the DNA fragments. 



With the increasing importance of biophysics there is a 
greater overlap between the subjects of the physics of bio- 
logical molecules and colloidal physics, which deals with 
small particles in suspension. Many of the techniques 
that have been developed in colloidal physics are directly 
applicable to biological molecules, such as proteins, and 
cells. DNA in particular is a very interesting biomolecule 
which exhibits a wide range of behaviors, due to its in- 
teractions with proteins and enzymes, but also due to 
its physical characteristics. A particular feature of DNA 
that is of interest to colloidal scientists is its high lin- 
ear charge density [l|, [2| . The interaction energy between 
DNA molecules and the structure of suspensions of DNA 
fragments have been shown to be strongly affected by 
this charge density. Consequently, other effects associ- 
ated with a strong charge, such as electrolyte screening, 
electrolyte friction, charge condensation, charge inver- 
sion and like-charge attraction in the presence of mul- 
tivalent salt ions, have also been predicted and investi- 
gated @, 0, 0, Q • Several investigations have considered 
the effect on the structure due to the renormalization of 
charge @,0|- Little work, however, has been performed 
on the interesting consequences of a high charge on the 
dynamics of these molecules. 

Highly charged colloids present an entirely different 
paradigm of particle interactions to the classical hard- 
sphere model. Well-known features of a low-salt suspen- 
sion of charged colloids are the low osmotic compressibil- 
ity of the suspension, and that the mean distance between 
the particles scales with the inverse of the cube root of 
the colloid concentration over a wide range of concentra- 
tions Q . This scaling corresponding to the formation of 
a correlation hole around the particles, leads to the dom- 
ination of the far-field over the near-field hydrodynamic 
interactions (HI), and therefore to an altered dynamic 
behavior which is unlike that of hard-sphere colloids. 

Wilk et al. [H have measured the long-time transla- 
tional self-diffusion coefficient, D L = (d^ + 2d£\ /3, 
of isotropic dispersions of 20 base pair oligomers of DNA 
by fluorescence correlation spectroscopy for various salt 
concentrations. DNA is a suitable molecule for study- 
ing effects resulting from electrostatic and hydrodynamic 



coupling due to its large linear charge density of approx- 
imately — 2e/3.4 A. The 20-mer DNA are almost per- 
fectly monodisperse, rigid cylindrical rods with length 
L — 6.8 nm and diameter d = 2 nm (aspect ratio of 3.4), 
a bare valency of Z = —42 in neutral or basic pH so- 
lutions and the translational free diffusion coefficient of 
D° = 1.07 x 10~ 6 cm 2 s _1 . The long-time coefficient, 
D L , was found to have an unexpected non-monotonic 
concentration dependence in low-salt conditions. 

In this letter we describe a versatile theoretical scheme 
that we have developed for the calculation of D L in 
colloidal systems 0, 0, [l(| • This scheme includes the 
long-range far-field part of the HI between the particles, 
which dominates in low-salt suspensions. We will show 
that in combination with colloid charge renormalization, 
this scheme successfully describes the non-monotonic de- 
pendence of D L ((f) on the macroion volume fraction <j>. 
This effect is of general importance and may be ob- 
served in any dispersion of colloids or biomolecules where 
long-range repulsive interactions are prevalent. The non- 
monotonicity in ip is unusual since it requires a delicate 
interplay of HI and electrosteric repulsions over a suffi- 
ciently broad concentration range. In contrast, a non- 
monotonic ^-dependence is not uncommon for transport 
properties such as the primary electroviscous coefficient, 
p(4>), associated with the suspension viscosity 77 [llj . 
and the collective diffusion coefficient D c 0]. Moreover, 
the electrophoretic mobility fi [l2j as well as 1/rj and 
D L (<j) = 0) exhibit a minimum as a function of the elec- 
trolyte concentration. The minimum in D L (0) is found 
for globular macroions [9|, LL3J , and also for semi-flexible 
charged polymers as shown recently in experiments and 
simulations In all cases considered, HI play a deci- 
sive role, e.g., the maximum of D c ((f>) arises from a bal- 
ance of the slowing HI and the speed-up of density relax- 
ations caused by the electrosteric repulsion. The maxi- 
mum in p(<p) at intermediate salinity arises from a compe- 
tition between the velocity gradient inside the macroion 
double layer that grows with 0, and the shrinking double 
layer distortion . Ignoring HI can lead to nonphysical 
results such as the failure to predict the maximum in the 
electrophoretic mobility of a short polyelectrolyte chain 
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as a function of the monomer number 15, 1 



Due to the increasing power of computers, substan- 
tial progress has been made in the simulation of trans- 
port properties in non-dilute charged colloidal disper- 
sions. Formerly, simulations that include the electroki- 
netic effect of electrolyte ions have been difficult due to 
the large asymmetry between these components, both in 
size and charge. These simulations have focused lar gely 
on the challenging problem of electrophoresis 12, 12, ljf . 
frequently used for particle characterization but for which 
a complete theory in dense systems, in which there is 
strong overlap of the electrical double layers, is still lack- 
ing. Our theoretical scheme is therefore also presented 
as a significant step in developing a versatile statistical 
mechanical tool to describe the electrokinetic transport 
in non-dilute suspensions. 

Our scheme is based on an exact memory equation for 
the self-intermediate scattering function of colloids un- 
dergoing overdamped Brownian motion. The irreducible 
memory function in this equation is approximated using 
the idealized mode-coupling theory for Brownian fluids, 
with the important extension of including the HI between 
all ionic species. Instead of solving the mode-coupling 
equations fully self-consistently, which would be a very 
challenging numerical task in the presence of HI, we use 
a simplified solution scheme which retains analytical sim- 
plicity and yields only small differences in the numerical 
results [lit . 

According to our scheme, the long-time coefficient 
D L of a non-dilute suspension is given by the Stokes- 
Einstein-like relation 0] 
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which includes, in addition to the colloid-solvent fric- 
tion C° = 67177 fl, where rf is the solvent viscosity and 
a the colloid radius, a colloid friction (CF) and an elec- 
trolyte friction (EF). The colloid friction arises from the 
microion-averaged electrosteric and HI between the col- 
loids, and is present even when the microion degrees of 
freedom are ignored. It is given in our scheme by [3j 
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where n is the colloid number density, h{k) is the total 
correlation function of the colloids, h d (k) is the distinct 
hydrodynamic function of the colloids, and D° is the free- 
diffusion constant. This expression for AC CF only re- 
quires the colloid static structure factor S(k) = l+nh(k). 
Eq. ([2]) is the zeroth order term in the expansion of the 
total long-time friction coefficient, AC = AC CF + A( EF , 
in terms of the colloid- microion mobility ratio D°/D®. 
Eq. H]) states that the extra friction due to the fast ki- 
netics of the mobile salt ions, A( EF , is given by the dif- 
ference of the total friction coefficient and A£ CF Q . The 



EF is due to the non-instantaneous relaxation of the mi- 
croionic atmosphere. We have derived a simple expres- 
sion for the long-time EF contribution, valid for the case 
when the free diffusion coefficients of the various salt ion 
species, D®, are much greater than the colloid free diffu- 
sion coefficient. This expression is 
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where the sum goes over all microion species of number 
density m, and where h C i(k) and hf^k) are the partial 
total correlation and partial distinct hydrodynamic func- 
tions between the microions and the colloids. The EF 
contribution to D L is significant in very dilute systems 
but is negligible when the mobility difference between the 
colloidal spheres and the microions is large, and when 
<j> is increased. The second finding is attributed to the 
enhanced homogenization of the electrolyte background 
with increasing <f> Q. 

Since the dynamic effect of the microions is small, we 
can simplify the problem by considering the colloids as 
interacting with an effective pair potential. In the result- 
ing one-component model of weakly charged colloids, the 
effective interaction between colloids of radius a and bare 
valency Zb aTe , follows the repulsive part of the DLVO po- 
tential si, 
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where Lb = e 2 / (AireksT) is the Bjerrum length in a 
solvent of dielectric constant e and K is the inverse De- 
bye screening length determined by the concentration of 
added salt ions and monovalent counterions. In a 1:1 
electrolyte solution, n 2 = AitLb (2n s + n\Zb alc \). For 
strongly charged colloids where Ls|Zb a rc|/a > 1, this 
potential is still suitable but only with the charge and 
screening parameter replaced by an effective charge, Z e g, 
and screening parameter, K c ff, due to the condensation 
of counterions near the colloid surfaces. There exist sev- 
eral schemes for the calculation of these effective quanti- 
ties, and these have recently been of considerable interest. 
Those mostly used are the cell model approximation of 
Alexander et al. [21 1, a nd the renormalized jellium ap- 
proximation (RJA) 22J. 

The RJA for the effective macroion charge in a closed 
suspension with a fixed salt concentration, as opposed 
to a system in contact with a reservoir, involves numer- 
ically solving the Poisson equation for the mean electric 
potential ip(r) of a single colloidal sphere, surrounded 
by a Boltzmann-distributed microion cloud and a uni- 
form negatively-charged background of charge density, 



3 



nZ e g, that describes the jellium representing the other 
macroions. The resulting equation for the reduced po- 
tential, y(r) = f3eip(r), is 

V 2 y(r) = AttLb 2n s sinh(y(r)) + nZ c g(e v(r) - 1) (3) 

with the boundary conditions y(oo) = 0, y'(oo) = 
0, and y'(a) = —LbZ^c/cl 2 . The numerical solu- 
tion of Eq. ([3]) is asymptotically matched to the solu- 
tion of the linearized equation, V 2 yn n (r) = K% s yn n (r), 
where n c g is the effective screening parameter k 2 s = 
AitLb {2n s + n\Z c g\). The effective charge comes from 
the linearized solution at the inner boundary, y[ in (a) = 
— LsZ c g/a 2 . Since the effective charge also appears in 
Eq. ([3]) this solution procedure must be iterated until 
self-consistency in Z c g has been established [jjj ]. 

With the so-determined interaction potential, S(k) is 
calculated using standard integral equation theories. For 
this study, we use the rescaled mean spherical approxi- 
mation (RMSA) [23| , known to be a reasonably accurate 
theory of the structure of charged colloids. 

Eq. ((2]) has been developed for spherically symmet- 
ric particles and is therefore not directly applicable to 
rodlike particles in high-salt suspensions. On the other 
hand, in the low-salt suspensions, as measured by Wilk 
et al., the microstructure is mostly determined by the 
long-range electrostatic monopole term which is spher- 
ically symmetric. Therefore, we treat the DNA frag- 
ments as spheres with an effective radius, a c g. This ra- 
dius appears in the solution of the effective charge and 
thus in the calculated structure factor. The largest ef- 
fect of a e g however, is to determine the scaling used 
to map the volume fraction dependent calculations onto 
the weight concentration dependent measurements by 
c(g/L) = 30M to /(4 7 ra3 ff lO 3 ^), where M w = 13022 g 
mol -1 is the molecular weight of the DNA fragments. 
The effective radius resulting from the Stokes-Einstein 
relation applied to the measured diffusion coefficient at 
infinite dilution is a c g = 2.0 nm. However, since the 
excluded volume interaction of the effective spheres is in- 
fluential at high salinity only, from comparing our results 
with the high-salt measurements of Wilk et at, we have 
determined a c g = 3.4 nm = L/2 to provide the best over- 
all fit, independent of the salt concentration. An effec- 
tive radius of half the molecular length has some parallels 
with the excluded volume calculations of rod-like parti- 
cles in the isotropic state. Since the bare charge is given 
by the number of ionizable groups on the DNA molecule, 
a c g = L/2 is the only adjusted parameter. 

The effective charge for our model of the DNA frag- 
ments, calculated via the RJA in a closed system as in 
the experiment, is presented in Fig. [T] For low-salt sys- 
tems, Z e g shows a non-monotonic dependence on (f>. At 
very low </), there is also a non-monotonic dependence of 
Z e g on the salt concentration (c.f., Fig. 1), with the limit 
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FIG. 1: Effective charge of 20-mer DNA (Z halc = -42), de- 
termined via the RJA, as a function of DNA and salt concen- 
trations. 



non-monotonicity of Z c g is also seen in its expansion in 
terms of the bare charge Zbare- If ^bare is sufficiently 
large, Z e g becomes independent of the bare charge. This 
is the so-called saturated effective charge, Z^g- In our 
system, \ZhaxeL b / des\ ~ 9 which is well into the non- 
linear regime, but Z e g is less than 66% of Z^g, so that 
saturation is not yet reached. Even though the K e g de- 
pends on Z e g, it shows no non-monotonic (^-dependence. 

The measured D L of the 20 -mer DNA and the com- 
parison with the results of our spherical model are shown 
in Fig. [21 In Fig. [5^, we include the results of our 
scheme when far-field HI are included (solid lines) or ig- 
nored (dashed lines), with the values Z e g and K e g used 
in both cases. The non-monotonic dependence of D L 
on <j> in the experiments is also seen in the theoreti- 
cal predictions when HI between the particles are in- 
cluded. If the HI are ignored, D L shows a monotonic 
(^-dependence, for all considered salt concentrations. In 
Fig. Wp, a comparison is made between the theoretical 
results for D L , including the effects of HI, where in- 
stead of a ^-dependent Z c g (solid lines again) we use 
a fixed charge of Z = —42 (dash-dotted lines). For both 
sets of results a non-monotonic ^-dependence is seen. 
The non-monotonicity, however, is much stronger when 
the non-constant Z e g is used. This suggests that the 
non-monotonic (^-dependence of the experimental D L re- 
sults from a simultaneous interplay between the hydrody- 
namic enhancement caused by far-field HI, and the non- 
monotonic Z c g in low-salt systems. The strong decline in 
D L at large <j) seen in Fig. [2] is due to electrosteric caging 
which becomes stronger with increasing (f>. The single- 
macroion EF effect described in Booth's theory 13( can 
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when n — > and n s 



0. This zero- 



not explain the non-monotonicity of D since it is sig- 
nificant only for <p < 10 -4 and n s w 0.01 M which 
is a salt concentration much larger than those where the 
DNA-D L behaves non-monotonically. 

The hydrodynamic enhancement of the diffusion of 
particles repelling each other over long distances results 



4 




tions. 

This work and the printing of the article was under 
appropriation of funds from the Deutsche Forschungsge- 
meinschaft (SFB-TR6, TP B2). We acknowledge A. Ban- 
chio, J. Gapinski and A. Patkowski for discussions. 



10 1 

DNA^ 



c~..(g/i) 



FIG. 2: Self-diffusion coefficient D L of the DNA model vs 
weight concentration for salt concentrations as indicated. 
Symbols are experimental results taken from [l[. a) Compar- 
ison with theoretical results using Z g. Solid lines are results 
with HI, dashed lines are without HI. b) Comparison with 
theoretical results with HI. Solid lines are again results using 
Z e s, dashed lines are results with fixed charge Z — —42. 



from the fact that the dominant far-field HI advect neigh- 
boring particles that may otherwise have hindered the 
motion of the considered one. Near-field HI, on the 
other hand, have the opposite effect of slowing the dif- 
fusive motion. Hydrodynamic enhancement of D L has 
been seen before in suspensions of charged and magneti- 



and in simu- 



cally interacting colloidal particles 1^, 
lations with HI of charged nano-sized polyions [25| , but 
without a visible non-monotonic 0-dependence. A non- 
monotonic D L ((f>) was found in simulations of salt-free 
polyelectrolyte solutions j26|. In these simulations, how- 
ever, HI have been neglected and the values for D L /D° 
are smaller than 0.1, which is the value where the freezing 
transition of charged spheres and rods occurs [27 1. 

In summary, we have shown that the non-monotonic 
concentration dependence of D L in low-salt suspensions 
of DNA fragments can be understood by the influence 
of far-field HI and charge renormalization. According to 
our scheme, the non-monotonicity of D L is a generic ef- 
fect for any low-salt suspension of strongly charged small 
colloids or bio-molecules. We have obtained this result 
using a simplified mode-coupling scheme. This scheme 
is a marked improvement on previous methods, since 
its many-component version includes the far-field HI be- 
tween all ionic species. Its analytic simplicity allows the 
study of electrokinetic phenomena such as the EF effect 
on self- and collective diffusion for non-zero concentra- 
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